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Abstract 

We construct a spectral representation for the baron to meson transition distribution am- 
plitudes (TDAs), i.e. matrix elements involving three quark correlators which arise in the 
description of baryon to meson transitions within the factorization approach to hard exclusive 
reactions. We generalize for these quantities the notion of double distributions introduced in the 
context of generalized parton distributions. We propose the generalization of A. Radyushkin's 
factorized Ansatz for the case of baryon to meson TDAs. Our construction opens the way to 
modeling of baryon to meson TDAs in their complete domain of definition and quantitative 
estimates of cross-sections for various hard exclusive reactions. 
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1. INTRODUCTION 



The concept of generalized parton distributions (GPDs) []l]-f§|, which in the simplest 
(leading twist) case are non-diagonal matrix elements of quark-antiquark or gluon-gluon 
non local operators on the light cone, has recently been extended |5], [J to baryon to meson 
(and baryon to photon) transition distribution amplitudes (TDAs), non diagonal matrix 
elements of three quark operators between two hadronic states of different baryon number 
(or between a baryon state and a photon). Nucleon to meson TDAs are conceptually much 
related to meson-nucleon generalized distribution amplitudes [|7|, |§ since they involve the 
same non-local operators P-pH- These objects are useful for the description of exclusive 
processes characterized by a baryonic exchange such as backward electroproduction of 
T~3| P~5]| or proton-antiproton hard exclusive annihilation processes fllBj . Nucleon 



mesons 



to meson TDAs are also considered to be a useful tool to quantify the pion cloud in 
baryons |17 . 



Up to now TDAs between the states of unequal baryon number lacked any suitable 
phenomenological parametrization in the whole domain of their definition, as for example 
in the framework of the quark model developed in | 18| . The complete parametrization 



should properly take into account the fundamental requirement of Lorentz covariance 
which is manifest as the polynomiality property of the Mellin moments in the relevant 
light-cone momentum fraction on the complete domain of their definition. For the case 
of the GPDs an elegant way to fulfill this requirement consists in employing the spectral 
representation. The corresponding spectral properties were established with the help of 
the alpha-representation techniques |19|, |2~0fl . Radyushkin's factorized Ansatz based on 
the double distribution representation for GPDs |2~l|-p4|1 became the basis for various 



successful phenomenological GPD models (see [^5| -[29|). 

In this paper we address the problem of construction of a spectral representation of 
baryon to meson transition distribution amplitudes. We introduce the notion of quadruple 
distributions and generalize Radyushkin's factorized Ansatz for this issue. This allows 
the modeling of baryon to meson TDAs in the complete domain of their definition and 
quantitative rate estimates in various hard exclusive reactions. 

Similarly as the nucleon to meson TDAs factorize in backward meson electroproduc- 



tion, nucleon to photon TDAs may factorize in backward virtual Compton scattering 



The main part of the analysis performed in our paper can be directly applied to the nu- 
cleon to photon TDAs. But the modelling of the quadruple distribution has to account 
for the anomalous nature of a photon. The studies of the anomalous photon structure 



functions pTJ and of the photon GPDs |32fl show that taking it into account is a nontrivial 



task which deserves separate studies. 

2. BASIC DEFINITIONS AND KINEMATICS 

Nucleon to meson transition distribution amplitudes also called in the literature as 
skewed DAs || and superskewed parton distributions || which extend the concept of usual 
generalized parton distributions arise e.g. in the description of meson electroproduction 
on the nucleon target [13-[R|. For defmiteness below we consider the case of nucleon to 



pion transition distribution amplitudes (nN TDAs for brevity) although our analysis is 
general enough to be applied to other baryon- meson and also baryon to photon TDAs. 
ttN TDAs arise in the description of backward pion electroproduction 

7*0?) + N{ Pl ) ^N\p 2 ) + Tt(p n ), (1) 

in the generalized Bjorken regime (— g 2 -large; q 2 /(2pi ■ q) kept fixed; — q 2 ^> —u). The 
factorization theorem was argued for the process ([]]) in || |(J (see Fig. [I]). The appropriate 



kinematics is described as follows [13) : 



P=^{Pi+Pn); A = p n -p l ; 

.9 A ■ n 
u = A ; f = , 2 

where u denotes the transfer momentum squared between the meson and the nucleon 
target and £ is the skewness parameter, n and p are the usual light-cone vectors occurring 
in the Sudakov decomposition of momenta (n 2 = p 2 = 0, n ■ p = 1). The light-cone 
decomposition of the particular vector v M is given by v M = v + p >1 + v + u^. 
The definition of ttN TDAs can be symbolically written as [§, |§[: 
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n dz i 

i=i 



e ixx(P-zi)+iX2(P-Z2)+ix3(P-Z3) 



x (7r(P + A/2)|e a6c ^ 1 (^i)V) 2 (^)Vi 3 (^)i^(P- A/2)}|^ =0 

~ £( 2 £ -x 1 -x 2 - x 3 )H jlhh (xi, x 2 , x 3 , f, u) . (3) 




FIG. 1: The factorization of the process 7* + P — > P' + ir. The lower blob is the 
pion-nucleon transition distribution amplitude, denotes the hard subprocess 
amplitude, DA is the nucleon distribution amplitude. 

Here 71,2.3 stand for spin-flavor indices and a, b, c are color indices. The decomposition 
of the Fourier transform (Q) of the matrix element of the three-local light-cone quark 
operator involves a set of independent spin-flavor structures multiplied by corresponding 
invariant functions: ttN TDAs. 

It is worth to mention that in order to preserve gauge invariance one has to insert the 
path-ordered gluonic exponentials [zf, zq] along the straight line connecting an arbitrary 
initial point zqu and a final one ZiW. 



(7T e 



a&c„/,a' 



i) a h {z X ) [21; Zo)a',a^ j2 (z 2 ) [z% Z )y b ^ (z 3 ) [z 3 , Z ]c'e \N) . (4) 



Throughout this paper we adopt the light-cone gauge A + = 0, so that the gauge link is 
equal to unity. Thus we do not show it explicitly in the definition (|3|). 

For the case of proton to tt° transition the decomposition of (§) over the independent 
spinor structures at the leading twist involves 8 independent terms. It reads 1 ||15|| : 

4J r ((7T (> 7r )| e abc u a a (z 1 n)u b f} (z 2 n)d^(z 3 n) \P{p u si))J = <5( 2 £ -x x -x 2 - x 3 ) 
xi fN \ V f^ CU{N+)i + °(^ 7 5 C) a/3 (7 5 iV + ) 7 + Tr\a Pfl CU(YN + ), 

+M- 1 vf{i>C) a ^N + ) 1 + M- x Af °(^7 5 C , ) a/3 (7 5 AriV + ) 7 + M^rf (a^C)^^), 
+M- 1 Tr\a PIX CU(a^N + ), + M- 2 Tf\a p ^ T C) a ^N + ), 1 . (5) 



1 We make use of the notation = (P ■ n) 3 J JlLi W e ix ^ p - z ^ +iX2( ~ p - z ^ +ix ^ p - z ^ (•) 
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Here j> is the usual Dirac slash notation {j> = p^j^), = |[7 M ,7^] with cr PM = p u cr u ^, 
C is the charge conjugation matrix and iV + is the large component of the nucleon spinor 
(N = {ftp 1 + j)ji)N = N~ + N + with iV+ ~ y/pf and N~ ~ y/l/pf)- M stands for the 
nucleon mass, f n is the pion decay constant (/„- = 131 MeV) and Jn is a constant with 
the dimension of energy squared. All the 8 p — > tt° TDAs V{, A{ and Ti are dimensionless. 

In this paper we concentrate on the dependence of the invariant functions Vi, Ai, Ti 
multiplying the independent spin-flavor structures in ([5]) on the longitudinal momentum 
fractions xi, X2, x 2 and skewness parameter £. Let us stress that our subsequent analysis 
is completely general: all invariant functions can be treated at the same footing. For 
simplicity in what follows we employ the same notation for all the invariant functions 

H(xt, x 2 , x 3 , f , t) = {V, Ai, Ti} (xi, x 2 , x 3 , £, u) . (6) 

A basic feature of model building cleverness is to fulfill fundamental requirements of 
field theory, such as general Lorentz covariance. In particular this requirement leads to the 
so-called polynomiality property of the Mellin moments in light-cone momentum fractions 
X\, x 2 , £3 of 71 N TDAs: 



J dxidx 2 dx 3 5(2£ — x\ — x 2 — x^x™ 1 x^ 2 x% 3 H {x\, x 2 , x 3 , £, u 
dz{ ) \ dz^ ) \ dz 3 



»'T^J ^T- ^ ) [( 7 r(P + A/2)|^(^i)V(^)V(^)|iV(P-A/2))] 



Zi=0 

(7) 



Indeed the X\, x 2 , x 3 - Mellin moments of irN TDA are the form factors of the local twist-3 
three quark operators between nucleon and pion states. This leads to the appearance of 
polynomials in £ at the right hand side of (^) 2 . 



2 Naive counting gives n\ +77,2 + 71.3 for the order of this polynomial. However, the problem of determina- 
tion of the highest possible power of £ in (0) still lacks some analysis. This is a rather important question 
since it would allow to make the conclusion on the necessity of adding of D-term like contributions J33| 
to the spectral representation of irN TDAs (see discussion in Sec. 
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3. SUPPORT PROPERTIES OF ttN TDAS 



3.1. ERBL-like and DGLAP-like domains for ttN TDAs 

In order to specify the support properties of ttN TDAs let us first consider the case 
of the GPDs (see Fig. 0.a). Let X\ and x<i be the fractions (defined with respect to 
average nucleon momentum P = 2i±£2) of the light-cone momentum carried by quark 
and antiquark inside nucleon (xi + 22 = 2£). In the so-called ERBL region both x\ and 
x 2 are positive. The variable x is usually defined as 



In the ERBL region x\, £2 £ [0, 2£] and thus x G [— £, £]. In the so-called DGLAP region 
either x\ is positive x\ G [2£, 1 + £] and 22 is negative 22 G [— 1 + £, 0] or vice versa (21 
is negative 21 G [—1 + £, 0] and X2 positive X2 G [2£, 1 + £]). These two DGLAP domains 
result in x G [£, 1] and x G [—1, — £] respectively. 



Xi Xi Xi x 2 x 3 




(a) (b) 

FIG. 2: Longitudinal momentum flow in the ERBL regime for GPDs (a) and nN 

TDAs (b). 

Now let us turn to the case of ttN TDAs. Let X\, 22 and 23 satisfying the constraint 
X1+X2+X3 = 2£, with £ > be the light-cone momentum fractions carried by three quarks. 
As usual the light-cone momentum fractions are defined with respect to the average hadron 
momentum P = P1 ~^ p ^ . The convenient way to depict the support properties of 7iN TDAs 
is to employ barycentric coordinates (Mandelstam plane). 

First of all we identify the analogous of the ERBL domain, in which three longitudinal 
momentum fraction carried by three quarks are positive. In the barycentric coordinates 
the ERBL-like region corresponds to the interior of the equilateral triangle with the height 
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2£ (see Figure |3[). It is natural to assume that the DGLAP-like domains are bounded by 
the lines 



Xi = — 1 + £; x 1 =0; Xi = l + £; 



x 2 = -1 + £ ; x 2 = ; x 2 = 1 + £ ; 



£ 3 = -! + £; x 3 = 0; a; 3 = l+£. 



(9) 



We are guided by the following requirements. 

• The complete domain of definition of ttN TDA should be symmetric in xi, x 2 , x 3 . 

• In the limiting case £ = 1 this domain should reduce to the ERBL-like domain 
on which the nucleon DA is defined. In the barycentric coordinates the domain of 
definition of the nucleon DA is equilateral triangle. 

• For any X{ set to zero we should recover the usual domain of definition of GPDs for 
the two remaining variables. 

Three small equilateral triangles correspond to DGLAP-like type I domains, where 
only one longitudinal momentum fractions is positive while two others are negative. Three 
trapezoid domains correspond to DGLAP-like type II, where two longitudinal momentum 
fractions are positive and one is negative. 

The support properties are invariant under the permutation of the longitudinal 
momentum fractions X{. In the limit £ — > 1 the support of ttN TDA is reduced to the 
ERBL-like domain (the equilateral triangle) (see Fig. [|) and coincide with that of the 
nucleon distribution amplitude (DA). In fact this is natural since £ = 1 corresponds to 
the soft pion limit in which irN TDA reduces to the corresponding nucleon DA ]T3] . 

In the limiting case £ — > the support of ttN TDA in the barycentric coordinates is 
given by the regular hexagon. 

3.2. Quark-diquark coordinates 

In order to describe irN TDA instead of x\, x 2 , £3 which satisfy 



xi + x 2 + x 3 = 2£ 



(10) 
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FIG. 3: Physical domains for ttN TDAs in the barycentric coordinates. 





FIG. 4: Physical domains for ttN TDAs in the barycentric coordinates. Two limiting 

cases: £ = (left) and £ = 1 (right). 



it is convenient to introduce the so-called quark-diquark coordinates. Let us stress that 
we do not imply any dynamical meaning to the notion of "diquark". There are three 
different possible choices depending on which quarks are supposed to form a "diquark 
system" : 



v 3 



x 2 




£3 




2 












2 




Xl 




x 2 



Wi 

w 3 



Xl 


- x 2 


- x 3 




2 




X 2 


- x 3 


— Xi 




2 




X3 


— X\ 


- x 2 



'IF 
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We suggest to introduce the notations £[, £ 2 and £ 3 for the fraction of the longitudinal 
momentum carried by the diquark: 

X2 + X 3 

2 

Xi + x 3 
2 

X\ + x 2 
2 

The variables x\, x 2 , x 3 are expressed through the new variables (O) as follows: 

x x = £ + w x ; x 2 = V! + ; x 3 = -v x + £[ ; 
x\ = -v 2 + £' 2 ; x 2 = £ + w 2 ; x 3 = v 2 + £' 2 ; 
x\ = v 3 + £ 3 ; x 2 = -v 3 + £ 3 ; x 3 = £ + w 3 . (13) 





- Wi 




2 




- W2 




2 






2 



(12) 



3.3. ERBL-like and DGLAP-like domains for 7riV TDA in quark-diquark coor- 
dinates 



Let us consider how the ERBL-like and DGLAP-like domains for ttN TDA look like in 
quark-diquark coordinates. Throughout the rest of this section we employ the particular 



choice of quark-diquark coordinates (]TTD : 



Xi-x 2 x 3 -xi-x 2 , , t,-w 3 
v = v 3 = ; w = w 3 = ; £ = £ 3 = — - — . (14) 

The generalization for the alternative cases is straightforward. 

The ERBL-like and DGLAP-like domains for nN TDA in quark-diquark coordinates 
([14]). are depicted on Figure [|. In these coordinates the ERBL-like region corresponds to 



the central isosceles triangular domain. It is bounded by the lines 

v = -? (zi = 0); v = i' (x 2 = 0); w = -£ (x 3 = 0) . (15) 

DGLAP-like type I regions correspond to three smaller isosceles triangular domains. Fi- 
nally, three trapezoid domains correspond to DGLAP-like type II region. 
For w G [—1, — £] DGLAP-like region is bounded by 

« = ! + (x l = l + and „ = -!-£ + £' (x 2 = 1 + £) ■ (16) 
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> 




FIG. 5: ERBL-like and DGLAP-like domains for tcN TDA in quark-diquark 
coordinates (0). Three lines: w = — £ and v = ±£' form the isosceles triangle which 
corresponds to ERBL-like region. Three smaller isosceles triangles correspond to 
DGLAP-like type I region. Three trapezoid domains correspond to DGLAP-like type II 

region. 

For w G [— £, £] DGLAP-like region is bounded by 

v = -i + z-? = + v = -g- 

« = v = l-t + t' (x 2 = -l+0- (17) 

For w G [£, 1] DGLAP-like region is bounded by 

v = -i + £-e (xi = -i + 0; « = 

i/ = £' ; u = i-£ + £' (x 2 = -i + 0- (18) 

One can easily check that for £ > the following inequalities are valid: 



(19) 



w < -f w > -£ 

; and < 

£' = £ occurs on the line u> = — £. Thus the whole domain of definition of ttN TDA in 
quark-diquark coordinates depicted on Figure can parameterized as follows: 

-1<w<1; -l + |£-f'|< v <l-|£_£'|. (20) 
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Let us briefly summarize our result. 



• w E [— 1; — £] with u G [£'; 1 — £' + £] or f G [—1 + £' — £; — £'] correspond to 
DGLAP-like type I domains. 

• w G [— 1; — £] and t> G [— 6'] corresponds to DGLAP-like type II domain. 

• w E [— £; £] with t> G [— £'] corresponds to ERBL-like domain. 

• WE [— £; £] with w G [£'; 1 — £ + £'] or t> G [— 1 + £ — — £'] correspond to 
DGLAP-like type II domain. 

• WE [£; 1] withw G [-£'; or u G f] correspond to DGLAP-like 
type II domain. 

• «)£[(; 1] with v E [£'; — £'] correspond to DGLAP-like type I domain. 
The Mellin moments of nN TDAs in xi, x%, £3 computed with the weight 

dxM 6( Xl+X2 + X3 -20 (21) 



are the quantities of major theoretical importance. In the quark-diquark coordinates (JT 
the corresponding integrals can be rewritten as 
/•!+£ 

/ dxidx2dx3 5(xi + X2 + £3 — 2£)x" 1 x\^ 'x^ 3 H (xi, X2, £3 = 2£ — x\ — X2) 
J-l-H 

1 r i-\z-e\ 

dw / dv{v + £') rai (-v + £')" 2 (w + O n3 H(w, v, . (22) 

•1 ■/ -i+[€— €'l 

4. SPECTRAL REPRESENTATION FOR vriV TDAS FROM THE SUPPORT 
PROPERTIES AND THE POLYNOMIALITY CONDITION 



The double distribution representation |21[-|24| was found to be an elegant way to 



incorporate both the polynomiality property of the Mellin moments and the support 
properties of GPDs. In the framework of this representation the GPD H is given as a one 
dimensional section of the double distribution (DD) f(a,(3): 

-1+1/31 



f l r l-\f3\ 

H(x, = / d P da8(x-p- a£)f(p, a) . (23) 

J — 1 J-i+\B\ 
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The spectral representation was originally recovered in the diagrammatical analysis 
employing the a-representation techniques |L9], pD[ . The spectral conditions \(3\ < 1 and 
I | < 1 — |/3 1 ensure the support property of GPD |x| < 1 for any |£| < 1. 

The polynomiality property of the Mellin moments in x which resides on the funda- 
mental field theoretic requirements (Lorentz covariance) is ensured by the fact that the 



x dependence of GPD in (|23|) is introduced solely through the integration path. In [[54 
it was pointed out that the relation between GPDs and DDs is the particular case of the 
Radon transform. It is worth to mention that the polynomiality property is well known 



in the framework of the Radon transform theory as the Cavalieri conditions [35 



Now we propose to invert the logic. From the pure mathematical point of view rep- 
resenting GPD as the Radon transform of a certain spectral density is the most natural 
way to ensure polynomiality property. Postulating the polynomiality property of GPD 
and the support property |x| < 1 one can put down the spectral representation (^) and 
unambiguously recover the spectral conditions < 1 and |a| < 1 — |/3|. Let us stress 
that this does not provide the alternative derivation of fl23|) since there is no way to show 
independently the support property \x\ < 1 of GPD. However we think that this line 
of argumentation justifies the use of the Radon transform ( P^D which is a rather general 
representation for a function satisfying the polynomiality condition with the restricted 
support in x as the building block for the spectral representation of multipartonic gen- 
eralizations of GPDs and in particular for tcN TDAs. In order to derive the form of the 
spectral representation for ttN TDA let us first consider the simple example of ordinary 
GPDs. 



4.1. Test ground: spectral representation for GPDs 

We are going to treat the example of usual GPDs in a slightly unusual way which we 
find more suitable for further generalization. Let us introduce the light-cone momentum 
fractions X\ and X2 of the average hadron momentum carried by the quark and antiquark 
respectively. The variables X\ and x 2 satisfy the condition x\ + x 2 = 2£. The support 
property in xi, x 2 is known to be given by 

-1 + £<Zi<1 + £; -l+£<x 2 <l + £. (24) 
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In order to write down the spectral representation for GPD we introduce two sets of 
spectral parameters Pi >2 , ai, 2 . The momentum fractions xi ;2 are supposed to have the 
following decomposition in terms of spectral parameters: 

xi = £ + Pi + aif ; x 2 = f + P2 + "2^ • (25) 

The condition xi + x 2 = 2£ can be taken into account by introducing two (^-functions 
8 (/3i + P 2 )8(ai + a 2 ). This allows us to write down the following spectral representation 
for GPD H(x u x 2 = 2£ - x u £) : 

H(xi, x 2 = 2£- xi, £) 

= / dflidai dp 2 da 2 8(xi — ( - Pi - ai()8(Pi + p 2 )8(ai + a 2 )F(Pi, p 2 ,ai,a 2 ) ■ 

(26) 

Here by Qi y2 we denote the usual domains in the parameter space: 

fii, 2 = {|/3i, 2 |<l; |«i, 2 | < l-|/3i, 2 |}; (27) 

and F(Pi, p 2 ,ai,a 2 ) is a certain quadruple distribution. 

The important advantage of the spectral representation (|26"D is that it is symmetric 
under the interchange of the longitudinal momentum fractions x\ and x 2 . Note that 
the spectral conditions (f27|) ensure the support properties ( |2~3| ) both in x\ and x 2 . The 
(rii, n 2 )-th Mellin moments in x±, x 2 of H(xi, x 2 = 2£ — xi, £) are polynomials of order 
ni + n 2 of £: 

/■i+e 

/ dxi / <ix 2 <5(2£ — x\ — x 2 )x T l 1 x^ 2 H (xi, x 2 = 2£ — x\, £) 
J-x+t J-iH 

= [ dPidai [ d/3 2 da 2 (£ + pi + aiO" 1 (£ + Pi + c^)™ 2 

x8(Pi + p 2 )8(ai + a 2 )F(Pi,p 2 ,ai,a 2 ) = P ni+n2 (0. (28) 

Now we are about to show that the spectral representation ( p6"D is equivalent to the 
usual Radyushkin's representation ( p3|) for GPDs in terms of double rather than quadruple 
distributions. For this issue we can lift the two superfluous integrations employing the two 
delta functions. In order to perform this in the astute way let us introduce the natural 
spectral variables a±, P±: 

«1 ± «2 a Pi ± P2 , on x 

«± = — g ' ^ = — 2 — ' ^ ' 
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It is also useful to to perform the related change of the variables in the (xi, x 2 ) space in 
the initial spectral representation (|26|). The corresponding natural variables are 



x- = Xl 2 X2 = a- + P-£ ; and x+ = Xl ~^ X2 = £ + a + + f3+£ . (30) 
Thus instead of using ( Effi) we switch to the natural variables and consider: 
H(xi, x 2 = 2£- xi, £) 

= - I df3\da\ I d(3 2 da 2 5(x- — /3_ — a-£)6(ft+)6(a+)F(l3i, fa, «i, 012) ■ 
2 Jq 1 Jn 2 

(31) 

The appropriate definition of the integration domain in (|31~1) after the change of the 
variables (|29|) require special attention. In particular, 

1 rl rl /-l-ljS-l 

dp x I dfa... = 2 d(3_ dp+... . (32) 

1 J-i J-i J-i+\$-\ 

Now since 1 — | > and hence — 1 + | fa | < the integral over j3 + can be easily lifted 
with no influence on the integration domain in a+, a_. The problem of definition of the 
integration domain in a + , a_ in principle is reduced to change of the variables in the 
integral 

/a rb 
da 1 / da 2 8{a\ + a 2 ) ... , (33) 
■a J-b 

where a = 1 — + b — 1 — — It is much simplified due to the fact that 
/3_|_ = and thus a — b = 1 — \(3-\. This gives 

pa pa ra— |a— 

dai / <ia:2<5(«i + ^2) ••• = 2 / da:_ / da + 5(a + )... . (34) 



J — a J —a J~ a+\ot-\ 

Now the integral over a + can be trivially performed with the help of 5-function again 
producing no additional restrictions for the integration domain in a_ and /?_. The final 
result reads 

H(xi, x 2 = 2£- xi, £) 

• 1 rl-\P-\ 

d(3_ / cfa_5(x_ -fa- a_0 2F (/3_, ar_, -ar_) . (35) 

— J-l+\P_\ y v ' 

/08-,a-) 

Certainly we just recovered the known Radyushkin's result for the double distribution 
representation of GPDs. 

Let us just make a short summary of the crucial points. 
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• We started from the spectral representation for H(xi, X2 = 2£—xi,£) as the function 
of the skewness parameter £ and of two longitudinal momentum fractions x±, x<i 
satisfying the condition x\ + x<i = 2£. The form of this spectral representation 
ensured the proper support properties in xi, X2 as well as the polynomiality property 
of the corresponding Mellin moments in x\ and x%- The spectral density was a 
certain quadruple rather than double distribution. 

• The constraint x\ + X2 = 2£ was taken into account by the introduction of two 
^-functions restricting the integration domain in the space of spectral variables. 

• The two superfluous integrations can be lifted with the help of two 5-functions. This 
requires the special attention to the integration domain in the space of spectral 
parameters. This problem can be most easily solved by switching to the set of 
natural variables 60th in the space of spectral parameters and x\, X2 space. 

• In our toy exercise lifting the two integrations does not lead to any special restric- 
tions on the remaining spectral parameters a_, /3_ and we just recover the usual 
Radyushkin's result for the double distribution representation of GPDs. 

• We find the spectral representation ( Effi) which is symmetric under the exchange of 
x\ and X2 suitable for the generalization to the multiparton case. The analysis of 
ttN TDAs with the help of the approach discussed above is presented in the next 
subsection. 

4.2. Spectral representation for ttN TDAs 

We are now about to apply the ideas described in the previous section to the case of 
ttN TDAs. Let us consider ttN TDA H(xi, X2, x% = 2£ — X\ — x 2 , £) as a function of 
light-cone momentum fractions X\, X2 and 23 carried by three quarks. The three light- 
cone momentum fractions satisfy the condition x\ + X2 + £3 = 2£. The support property 
in xi, X2, X3 is given by 

-l + f<xi <! + £; -1 + £<x 2 <1 + £; -1 + £<x 3 <1 + £. (36) 
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In order to write down the spectral representation for H(xi, X2, x 3 = 2£ — x% — X2, £) 
we introduce three sets of spectral parameters 0i^,3, 01,2,3- The momentum fractions £1,2,3 
are supposed to have the following decomposition in terms of spectral parameters: 



xi = £ + 0i + cti£ ; x 2 = £ + 02 + a 2 £ ; x 3 = £ + 3 + a 3 £ . 



In order to satisfy this constrain we require that 



0i + 02 + 03 — ; a 1 + a 2 + a 3 = -1 



This allows to write down the following spectral representation for ttN TDAs: 



(37) 



(38) 



H(xi, x 2 , x 3 = 2£-x 1 - x 2 , 
3 



6(xx -£-0i- a x £) 5(x 2 - £ - 2 - a 2 f) 



XS(0! +0 2 + 0s)5(ai +CX2 + CX3 + l)F(0i, 02, 03, Oil, «2, oi 3 ) . 



By fli, i = {1, 2, 3} we denote the usual domains in the parameter space: 



n l = {\0 i \<i; H<i-|AI}; 



(39) 



(40) 



and F(0i, 2 , 03, a x , 02,03) is now a sextuple distribution. The spectral conditions (|40| ) 
ensure the support properties (|36|). Obviously, the (m, ri2, n 3 )-t\i Mellin moment in 
(xi, X2, x 3 ) of 77 N TDA is a polynomial of order n x + n 2 + n 3 of £: 



' 3 

n / 

i=i 



5(xi + x 2 + x 3 - 2£)x r l 1 X2 2 x 3 3 H(x 1 , x 2 , x 3 = 2£-xi~ x 2 , £) 



P, 



n 1 +n 2 +n 3 



(0 



(41) 



In complete analogy with the previously considered case of usual GPDs in order to 
properly reduce the spectral representation in terms of sextuple distribution for ttN TDA 
to that in terms of quadruple distribution we need to perform two integrations in 

3 „ 

5(0i + 02 + 3 )5(a 1 + a 2 + a 3 + 1) ... (42) 



] [ / d0idoii 

i=l 



employing 5-functions and specify the integration limits in the remaining four integrals. 
This problem can be solved by introducing the appropriate natural variables. 
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Let us start with the integral 

-i r i r± 



J d?! J dft J , 



dft5(ft + ft + ft). (43) 



In order to visualize the integration domain (33) it is natural to employ the barycentric 



coordinates. In these coordinates the domain selected by the conditions |ft| < 1 (i 6 
{1, 2, 3}) and ft + ft + ft = is represented by a regular hexagon (confer Fig It is 
convenient to single out three domains inside this hexagon: 

D 1 : {ft > 0, ft < 0, ft < 0} U {ft < 0, ft > 0, ft > 0} ; 
^2 : {ft > 0, ft < 0, ft < 0} U {ft < 0, ft > 0, ft > 0} ; 
£3 : {ft > 0, ft < 0, ft < 0} U {ft < 0, ft > 0, ft > 0} . (44) 

Obviously 

/ tZft f dft / dft5(ft + ft + ft) = V / dftdftdft5(ft + ft + ft). (45) 

./-l ./-l ./-l i=1 J Di 



Now in order to get rid of one of three integrations in (f45|) we should switch to the 
natural coordinates. There are three possible choices of the natural coordinates in (f4~5|) . 
For the moment we are going to adopt the coordinates 

ft -ft ft ~ ft - ft (AR , 
P3 = — 2 ' a3 = 2 ' ^ ^ 

The constrained triple integral ([43"D can be then rewritten as 



1 rl- 2 



l£3 

do-s / 2 rfp 3 -. (47) 

-1 

In principle in a completely analogous way one may also employ the coordinates 



ft - ft ft - ft - # 



pi = — ^ — ; °"i 



3 



2 2 



ft -ft ft - ft - ft ... 

P2 = „ ; 0"2 = ~ (46 



yielding the result 



1 /■! |CT 



2 



d<7j / dpi.... (49) 

1 y-i+isi 
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Now let us address the problem of computation of the constrained triple integral over 

a,- in my. 



/a rb t*c 

dcti / da 2 / da 3 5(ai + a 2 + a?, + 1) 
a J — b J —c 



(50) 



where we introduced the notations 

a = 1-1/3x1; b=l-\p 2 \; c = 1 - \/3 3 \ ; (51) 



a > 0; a < 1, 6 > 0; 6 < 1, c > 0; c < 1). 
Introducing the natural coordinates 3 



u 3 = a 3 ; j/ 3 = (52) 

and employing the results of the Appendix |A] we conclude that for fa G D\ U D 2 U -D3 the 
constrained integral ([50]) can be rewritten as 

du; 3 / dv 3 .... (53) 

i+|fe| J-i+m+^p- 

Now let us put all together and write down the spectral representation for 7riV TDAs in 

terms of quadruple distributions. The important observation is that once we have chosen 

the variables 03, p 3 and w 3 , z/ 3 to perform the constrained integration in fa, fa, fa and 

«i, a 2 , «3 respectively the natural variables on which ttN TDAs depends are 

x?,-xi-x 2 x 1 -x 2 
= g ' Vs = 2 — ' ^ ' 

Expressing the fa and «j through cr 3 , p 3 , w 3 , z/ 3 the two delta functions in the definition 
fl3~9l) can be traded for 

5(X X - £ - fa - «x0 6{X2-Z-P 2 - «20lxi+x 2 +x 8 =2{ 

= 5(w 3 - cx 3 - w 3 £) 5(^3 - P3 - ^0 • (55) 

Note that at the level of delta functions we achieved the "factorization" of W3 and i> 3 
dependencies on the spectral parameters. 



There are two additional possible choices: toi = olx ; v\— 012 a3 and W2 = «2 ; ^2 = ° 3 7 ai . 
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Thus in the natural spectral parameters ( ^6[ ) , ([)2|) and quark-diquark coordinates (|5 
we recovered the form of the spectral representation of tcN TDAs in terms of quadruple 
distributions: 

H(w 3 , v 3 , 

df3 1 d(3 2 d(3 3 5{Pi + (3 2 + (3 3 ) / c?q;i / rfa 2 / da 3 5(ai + a 2 + a 3 + 1) 

'-1 ^-l+|/9 2 | J-l+\Ps\ 

5(x! a^)5(x 2 -{-fa- afcO Wi, Aj, ai, a 2 , a 3 ) 

»1 /-1-|P3-^|-|P3+^| /-|-|P3+¥|-^ 

dcx 3 / rfp3 / du> 3 / du 3 5(w 3 - o- 3 - 



x5(w 3 - p 3 - f 3 £) F 3 (o- 3 , p 3 , w 3 , z/ 3 ) , (56) 
where 

7-1 / \ 7-1/ ^3 ^3 1 + W 3 1+W 3 

F 3 (cr 3 , p 3 , w 3 , v 3 ) = F(p 3 - — , -p 3 - y, a 3 , f 3 - — , -v 3 — ,u 3 ) . (57) 

Employing three possible sets of natural spectral parameters one can write down three 
equivalent spectral representations in terms of three sets of quark-diquark coordinates Wi, 
Vi with i — 1, 2, 3 defined in (|TTD : 



1 r i-*P f i-\ Pi -^\-\ Pi +?i-\ f L-\ Pi+ ?l\-»< 



H(w h Vi, £) 

doi \ dpi I dui I dvi5(wi - ai - 0Ji£) 

x8(vi - pi - Uig) Fifa, pi, Ui, i>i) , (58) 

where F 3 (a 3 , p 3 , u 3 , z/ 3 ) is defined in fl5"T|) and 

Fi(ai, pi, wx, z/i) = F(eri, pi - y , -pi - y , Wi, z/i — , -v x - — ) ; 

, , _ . a 2 a 2 1 + uj 2 l + u 2 , 
F 2 (a 2 , p 2 , ui 2 , v 2 ) = F(-p 2 - y , o- 2 , p 2 - —, -v 2 — , u 2 , v 2 - — ) . 

(59) 

The spectral representation (^) for irN TDA in terms of quadruple distribution is the 
main result of our paper. However this form of the result is still not very convenient for 
practical applications. In the next section we demonstrate that the spectral representation 
(|58|) satisfies the support properties of ttN TDAs established in Sec. ||. We also derive 
the explicit expressions for ttN TDAs in the ERBL-like and DGLAP-like type I and II 
domains. 
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5. SUPPORT PROPERTIES OF irN TDAS AND THE SPECTRAL REPRESEN- 
TATION 



In order to make our formulas more compact in what follows we omit the indice i for 
the quark-diquark coordinates Wi and Vi, spectral parameters <7j, pi, Ui, Vi and the spectral 
densities Fj. Our subsequent analysis equally applies for all i — 1, 2, 3. 

It is extremely instructive to check that each contribution into ttN TDA in (|58|) satisfies 
the support properties which were established in Sec. || 



with £' defined in (|12|). In particular this allows to check that (N — n, n)-th (N > n > 0) 
Mellin moments of ttN TDA in (w, u) indeed satisfy the polynomiality property: 

dw dvw N - n v n H(w, v, = P N (g) , (61) 

J —l J-i+\e-e\ 

where Pn{Q is a polynomial of order N in £. 



Let us first consider the case ^ = 0. Employing the first delta function we get a — w 
for — 1 < w < 1 and otherwise. This obviously ensures the first condition (|60|) for £ = 0. 
Once the integral over a is performed the dependence on v is introduced through 



- 1< w < 1; 



(60) 



Case £ = 




(62) 



The result of this integral is non-zero only for 



w 



(63) 



- 1 + — < v < 1 
2 - ~ 

that is precisely the second condition flSTfr ) for £ = 



0. 



2 



Case < £ < 1 



Let us now show that the spectral representation flSBp possesses the desired support 
properties for arbitrary value of £ G (0; l] 4 . 



4 The final result for £ £ [—1; 0) is presented in the Appendix |b| 
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First of all it is easy to see that the first one of the two conditions ( j6Qj ) is respected. 
Indeed the w dependence in (ESI) is introduced through the expression 



da 



i-|p-f|-b+f| 



dui 8{w — a — cu£) ... . 



(64) 



en 



(65) 



-1 J-l+\a\ 

From the inequalities flftSJ), (|A11| ) and ( |A14|) it follows that 

Thus in ( |6"4j ) we are integrating only over some part of the familiar "GPD square" |p| < 
1 — \a\. This guarantees the vanishing of ttN TDA for \w\ > 1. One can in the usual 
way perform the integration over u introducing the additional ^-function to take into the 
account the unusual upper limit in the integral over u: 

w — a , 



6(1 -\p 



a 

r 



a 
V 



9{...). 



(66) 



For £ > we get 



H(w, v, = 
For w G (— oo; —1) : ; 
For w G [-1; -£] : 



1 f 



do" 



For u> G [-£; f] : 



a_ I w — a 

2 25~ 



U> — cr 



dp / di>5(v - p- v£)0(...) F(a,p 7 

1 I [ _ £T I | If — (7 r" 



l+{ 1 ^ 2 

For iu G [£; 1] : 



/ 2 1^ 2 I 2? w — a 

dp dv8(v - p-vg)6(...)F(<T,p, — - — . 

/ 1 I „ ct I I w — a r 



1 f l+« 



da 



dp 



1+ 1 



~2+\p- 



a_ I v: — (T 

2 I 2£ 



rfz/5(w - p - I/O 0(...) F(<r, p, W - a , i/) ; 



ct I I w — a 
2 I" 1 " 2£ 



For w G (1; oo) : . 



(67) 



Now we are about to perform the integration over v with the help of the last remaining 
5-function. The resulting domain of integration in a and p is defined by the inequalities 



-1 



en 



< p < 1 



\a\ 

Y 



(68) 



21 



i-\p-- 2 \-\p + - 2 \>— ; (70) 




as well as the integration limits in a depending on the value of w (see 

It can be shown that for £ > the two inequalities ( |69|) are equivalent to 

P < \ + ^| for v > -£'; p < °- + for i; < (71) 

together with 

p > -- + ^| for « > P > -| + ^| for « < £' . (72) 

Analogously the inequality (^) for £ > is equivalent to 

P<^ + | for v>|£'|; p>-^-| for v<-|£'|; 

a>^ for r-~ r ; a > ^| for < " ' . (73) 

The last step is to match the integration domain defined by the inequalities (0), (|7l|), 
(0) and ( |73"D with the explicit w-dependent limits of integration in a (|67D. There are 9 
possibilities: 

{we[-l; -£], iu e [-£; -£], w e [£; 1]} 

®{ vG (-oo;-|ri], wGhiam «G[|ri;oo)}. (74) 
Let us consider in details the case 

we[-l;-£]; ue[e';oa). (75) 

The integration domain in (a, p) plane is defined by the intersection of a domain specified 
by the inequalities (0), (^), ([§: 

c f — C a v + £' a £' , , IcH , , 



with the strip 



^|<a<^±|. (77) 
22 



The domain denned by the inequalities (ITS) and (|77| ) is presented on Fig. |6|. By the thick 
solid lines we show the borders of the domain defined by the first two inequalities (|7"6|). 
The thin solid line is the border of the domain defined by the inequality p < 1 — | . The 
dashed line is the border of the domain defined by the inequality ([75]): p < + j- The 
shaded area corresponds to the resulting domain of integration in ( 
and f < v < 1 - £' + £. 

P 



for — 1 < w < — £ 




cr 



FIG. 6: The domain of integration in (a, p) plane in eq. ( |67"D for — 1 < w < — £ and 
£' < v < 1 — + £ defined by the inequalities (|76| ) and (|77|). See explanations in the text. 



The abscissa of the apex of this triangular domain is 



a 



1 -e 2 ' 

One may check that for v = £' the abscissa of the apex coincides with the left 
of the strip (fT?)): 



2K-0 



^ 2 



(78) 
boundary 

(79) 



while for v — 1 — £' + ^ it coincides with the right boundary of the strip (|77[): 



2K-0 



1-f 



1-^ 



?0) 



For u > 1 — £' + ^ the apex of the triangular domain lies on the right of the 
and hence has empty intersection with it. This makes the double integral (|67|) 
v > 1 — £' + ^ and ensures the desired support property of H(w, v,£). 



strip Q77|) 
vanish for 
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The third inequality in ( |76|) does not further restrict the domain since the apex of the 
triangular domain belongs to the line p = ^ + j and the triangular domain lies to the 
right of this line for < £ < 1. The two inequalities ( |B5| ) also do not impose additional 
restriction for the domain. Indeed one may check that the p = ^ + j intersects with 
p = 1 - M at a = f±f . 

The eight remaining cases ([74]) can be considered according to this pattern in a com- 
pletely analogous way. One may check that the quadruple integral (|58[) for H(w,v,^) for 
£ > reduces to the following expressions: 

• For w and v outside the domain w G [—1; 1] and v G [—1 + |£ — £'|; 1 — |£ — £'|] the 
integral vanishes. 

• For w G [-1; -£] and v G [£'; 1 - £' + £] (DGLAP-like type I domain): 

/ !-£ / 2 i+£ w — a V — 

H(w,v,0 = - da ^pF^p,^—,—?). (81) 

• For w G [-1; -£] and v G [-£'; £'] (DGLAP-like type II domain): 

H{w,v,i) = - da dpF(a,p,^-^/—P). (82) 

• For to G [-1; -£] and v G [-1 + £' - £; -£'] (DGLAP-like type I domain): 

/ 2 i-i w — a v — o 

H(w, „,£=_/ da / dp F a, p, — — , —£■). (83) 

• For w G [-£; £] and v G [£'; 1 - £ + £'] (DGLAP-like type II domain): 

/ i+€ / 2 W — (J V — 

H(w, v, £ = - / dcr / dp F(<t, p, -— -, — ^ . (84) 

• we [-£; £] and u G [-£'; £'] (ERBL-like domain): 

/ i+{ / 2 i+s w — a v — o 

H{w,v,i) = - da d pF{a,p,—?-,—^). (85) 

• to G [-£; £] and u G [-1 + £ - £'; -£'] (DGLAP-like type II domain): 

-uj+£ ct _i_ u+£ / 

H{w, v^) = - da / dpF(a, p, — — , — ^) . (86) 
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• w £ [£; 1] and v £ [— 1 — £ + £']: the result coincides with fl34| ) as it certainly 
should be since this is the part of the same DGLAP type II domain. Note that 
this makes H(w, v, £) a smooth function for w = £ as it should be since this line 
(wi = £ -v=> Xi = 2£) does not correspond to any change of evolution properties of 
H(w, v, 0- 

• w £ [£; 1] and u £ [£'; -£'] (DGLAP-like type I domain): 

H(w,v,0 = - 2 da ipF ^p,^—^^—±). (87) 

• ty £ [£; 1] and t> £ [— 1 + £ — £']: the result again coincides with fl86|) since this 
is the part of the same DGLAP-like type II domain. 



6. RADYUSHKIN TYPE ANSATZ FOR ttN TDAS 

In this section we discuss what could be a possible approach for modelling of quadruple 
distributions F(a, p, u, v) occurring in the spectral representation ([58]). 

Employing the analogy with the case of usual GPDs one may assume that the profile of 
F(a, p, co, v) in (cr, p) space is determined by the shape of the function f(a, p) to which 
ttN TDA is reduced in the limit £ — > 0. For the moment we put aside the complicated and 
interesting problem of the rigorous physical meaning of this limit. It will be discussed 
elsewhere. Thus, we suggest to employ the following factorized Ansatz for quadruple 
distributions: 

F(a, p, u, v) = f(a, p) h(a, p, u, v) , (88) 

where h(a, p, u, v) is a profile function normalized according to: 

i-IP-f |-[P+f I /-i-lp+fh 1 ^ 

duo / dvh(a, p, u, v) — 1 . (89) 

-i+kl ^-i+IP-fl+H^ 

A possible model is to exploit further the analogy with the standard Radyushkin Ansatz 

for the double distributions p4] and to assume that the (u, v) profile of h(a, p, u, v) is 



determined by the shape of the asymptotic form of the nucleon distribution amplitude: 

15 

$ as (2/i, 2/2, ys) = -j ymy 3 ■ (90) 
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(91) 



The DA @g) is defined for y 1|2 , 3 G [0; 2] such that j/i + y 2 + y 3 = 2. 

In terms of quark-diquark variables u = 1 — y\ — y 2 and v = yi ~ V2 $ as reads: 

Note that 

^ d yi dy 2 dy 3 S(2 - 2/1 - 2/ 2 - J/ 3 )$ as (2/i, y 2 , y 3 ) = f ^ f '_ e ^"(w, *>) = 1 (92) 



$ as (u), P) is defined for 

— 1 < a) < 1 ; and 
while h(cr, p, oj, v) is defined for 



1 — Co 1 — to 

< V < 



(93) 



-l + k|<a;<l-|p-||-|p+|| 



1 -u 



2 +lp 2 



< v < 



1-co 



\p 



2 2' 

Thus it makes sense to employ the following substitution of the variables: 



(94) 



v 



"+§( 


P-fl 


+ |P+f 




H). 


i-H 
(i-«) 


P-fl 
2i/ - 


+ 
P 


|P+f 

-ih 


l + H)' 
> + fl 


2 






p-fl 




p + fl 



(95) 



This results in the following expression for the profile function /i(er, p, u, v) 
h(a, p, uj, u) 



15(l + 2i/-w-2 


^"f 1) (!" 


2i/ 


— — 


2 


p+tDa- 


(7 




16 (i 




P-f 


+ 


P+f 


+ 


H)) 5 



.(96) 



One may check that the profile function (PB) satisfies the normalization condition (|89|). 
It is extremely interesting to note that in terms of the initial spectral parameters a%, a 2 , 
a 3 and (3%, (3 2 , (3 3 satisfying a.\ + a 2 + a 3 = — 1 and 0i + f3 2 + {3 3 = the profile function 
(PB]) can be rewritten in the very symmetric form: 

15 n-=i(i+«i-iAi) 



h(8 h 8 2 ,8 3 :ai,a 2 ,a 3 ). El(9l =o , , , 

Ei«i=-i 4 (1 - i ( |/5i | + |/3 2 | + |A?|)) 5 



(97) 
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The inverse transformation 
1 



reads 



u 1- 



OJ 



v = v\ \ 



a 




a 




a 




a 




P ~2 


+ 


" + 2 




P ~2 


+ 


" + 2 


-M) 



0" 






o 






(7 






a 




P "2 


+ 




' 2 






P "2 






" 2 


) 



(9? 



This allows to easily compute the integrals occurring in the calculation of (N — n, n)-th 
Mellin moments (N > n > 0) in (w, v) of ttN TDAs: 



i-|p-f|-|p+f| /-l-lp+fh^ 



<iu; 



i+H 



i+|p-f|+^ 



rfz/w n u n h(a } p, oo, v) . 



(99) 



In principle one may also think of a more intricate profile function. In fact any partic- 
ular function $(cD, v) normalized according to 



duo 



dv$(£j, u) — 1 



1— s 
2 



(100) 



will define some profile function h(a, p, u, v) after the substitution (p5|) 5 . E.g. taking 



$(w, z>) ~ (a) + \) bx {v + - L ^) b2 (-z/ + would lead to the natural generalization of 

the b parameter dependent Radyushkin profile familiar for usual GPDs. 

It is interesting also to consider the most simple possible profile with no distortion in 
(cj, v) directions: 



$(u), v) = S(Cj)S(u). 



(101) 



Contrary to the case of usual GPDs for which the counterpart of the profile (|101| ) leads 
to ^-independent Ansatz the resulting ttN TDA preserves the minimal necessary £ de- 
pendence. Indeed Cj = and v = does not imply oj = and v — and hence the 
^-dependence introduced through two 5-functions in ( |58f) is preserved and generates the 
proper ^-dependent domain of definition for the resulting ttN TDA (|20"D . Unfortunately 



the model with the profile (|101|) turns out to be pathological since it leads to ttN TDAs 
which are not continuous at the cross-over lines v = ±£' and w = — £ separating ERBL- 
like and DGLAP-like type I, II domains. This makes impossible the calculation of the 
amplitude of the hard exclusive process in question given by convolution of ttN TDA 



5 However, one has to make certain assumptions on the endpoint behavior of the function /(cr, p) to 
which irN TDA is reduced in the limit £—5-0. 



27 



with the appropriate hard part (see \TE\). Indeed the imaginary part of the corresponding 



amplitude is given by the values of irN TDA at the cross-over lines v = ±£' and w = — £. 

For the moment as a toy model we are going to employ the factorized Ansatz Q8"8]) 
with the profile function ([)6|). It is a good point now to discuss a possible model for the 
function f(cr,p) that is the second ingredient of the factorized Ansatz (pSf). In the limit 
£ — > nN TDA reduces to this function: 

H(w,v^ = 0) = f(w,v). (102) 



The requirements of convergence of integrals (|81|)-(p7D for irN TDA impose some 
restriction on the behavior of the function f(a, p) on the border of its domain of definition. 
It turns out that f(a, p) should vanish at least as a certain power of the relevant variables 
at the borders of its domain of definition. Thus for the function f(a, p) we suggest the 
following simple form: 

/!>, p) = < <J < 1) 9(-l + M < p < 1 - M) 



^t 1 -)^- 1 ) -Tjr +1) -tJ- (io3) 

In terms of the initial spectral parameters /3j satisfying £\ = (|103|) can be rewritten 
as: 

40 3 

/(A, A, /3 3 )l Eift=0 = ^ 11^1^1 ^ - (104) 

i=i 

The function /(a, p) vanishes on the border of its domain of definition and is normalized 
according to 



i 



da / , f(a, p) = l. (105) 



Let us stress that we employ the normalization ( p.05|) only for our toy model. Advanced 
modelling of ttN TDAs aiming the quantitative description of the physical observables 
would certainly require more complicated form of f(a,p). 

The normalization for the nucleon to pion TDAs can be derived either from the soft 
pion limit or from the lattice calculations of several first Mellin moments of ttN TDAs 
or from the comparison with the results of |[L8| . On the other hand it can be computed 
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considering the light baryon exchange contributions into the Mellin moments of 7rN TDAs 
using the phenomenological values say of g n NN and g n n a couplings. The normalization 
can also in principle be established directly form the experimental measurements of the 
cross-section once the scaling behavior would be found reasonable. 

On Figure we show the results of the calculation of the contribution H(u>3, V3, £) = 
H(w,v,£) as a function of w and v for different values of £ computed with the help of the 
factorized Ansatz ( j88|) with the profile (fjop and f(a, p) given by the toy model (|103|) . 

Note that for £ = 1 the TDA H(w,v,£) does not vanish at the corners of its domain 
of definition. This is potentially dangerous since this may lead to the break up of the 
factorization property of the hard exclusive process in question. Fortunately this problem 
is an artefact of our oversimplified toy model (|104j) for the forward limit of tcN TDA. It 
was checked that taking f(a,p) that vanishes quadratically at the borders of the domain 
of definition 

3 

m, /3 3 )i Eift=0 = ^ f[em\ < m-m\ if)i=0 (ice) 

i=i 

leads to a vanishing ttN TDA at the corners of its domain of definition for £ = 1. 

On Figure [8] we show ttN TDA H(xi,X2,Xa, £) for £ = 0.5 as a function of three 
dependent light-cone momentum fractions xi, £2 an d £3 {x\ + x<i + X3 = 2£) in the 
barycentric coordinates. By thick solid lines we show the continuation of the edges of the 
equilateral triangle which form the ERBL-like domain cf. Fig. [3]. 
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FIG. 7: The contribution into ttN TDA H(ws,vs,^) = H(w,v,£) as a function of u> 
and v for different values of £ computed using the factorized Ansatz fl88|) with the profile 
function ( p^) and /(a, p) given by ( |103| ). 
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H(x\,x 2 ,xj, £=0.5) 




FIG. 8: ttN TDA H{x\, X2, X3, £) as a function of x%, X2 and X3 {x\ + £2 + ^3 = 2£) for 
£ = 0.5 in barycentric coordinates. By thick solid lines we show the continuation of the 
edges of the equilateral triangle that border the ERBL-like domain cf. Fig. |3|. 

7. CONCLUSIONS 

The non-perturbative part of hard processes involving hadrons is encoded in various 
universal partonic distributions (parton distribution functions, fragmentation functions, 
distribution amplitudes and their generalizations). Waiting for a complete understanding 
of the dynamics of quark and gluon confinement in hadrons, one should model these 
distributions in agreement with general requirements of the underlying field theory such 
as Lorentz invariance and causality. Spectral representation of hadronic matrix elements 
offers an elegant way to address this program. The double distribution representation for 
GPDs became the basis for various successful phenomenological GPD models. 

In this paper we introduced the notion of quadruple distributions and constructed 
the spectral representation for the transition distribution amplitudes involving three par- 
ton correlators which arise in the description of baryon to meson transitions. We also 
generalized Radyushkin's factorized Ansatz for the case of quadruple distributions and 
provided an explicit expression for the corresponding profile function. Analogously to the 
case of GPDs the shape of the corresponding profile function is supposed to be fixed by 
the asymptotic form of the nucleon distribution amplitude. Our model also requires the 
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knowledge of nucleon to meson TDAs in the forward limit as input quantities. Contrarily 
to the GPD case, the nucleon to meson TDAs suffer from the fact that there is no illu- 
minating forward limit. This problem requires further investigation. For a moment, we 
suggest to employ a simple shape of nucleon to meson TDAs in the forward limit assuming 
that they are fixed by their behavior at the borders of their domain of definition. Our 
construction opens the way to quantitative modeling of baryon-meson and baryon-photon 
TDAs in their complete domain of definition. 

Let us emphasize that for the moment we have not included any D-term like contribu- 
tions to the spectral representation of the nucleon to meson TDAs in terms of quadruple 



distributions. Indeed the results of |36| and of Chapter 3.8 of [^j give us confidence that 
the eventual D-term like contributions to TDAs can be included by means of comple- 
menting the spectral density in (|39"D with additional terms proportional to powers of £. 
The subsequent analysis can be performed according to the same pattern. 

Let us also point out that our method can be generalized for the case of 4-quark 
correlators important for the description of higher twist contributions. 
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A. A USEFUL CONSTRAINED INTEGRAL 

Let us consider the constrained triple integral 

/a rb pc 

dai / da 2 / da 3 8{a,\ + a 2 + a 3 + 1) a 2 , a 3 ) , (Al) 
-a J —b J —c 

where a > ; a < 1, b > ; b < 1, c > ; c < 1. We introduce the natural coordinates u 3 
and v 3 : 

— 1 — U 3 — 1 — UI3 

"i = ^3 H ^ ; ol 2 = -v 3 H o ; a 3 = u 3 . (A2) 
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In the natural coordinates u 3 and v 3 the integration in ( |A1| ) is over the intersection of 
three stripes: 



— c < u> 3 < c ; 

. 1 + w 3 . 1 + w 3 
-a H z — - Us - a H 2 — ' 

-b- l -±^<v.<b- l A^. (A3) 



One may check that for a > b the integral (A~I) can be rewritten as 



/— 1— a+6 /-oH 2~ ^ 

ciw 3 6»(a; 3 + c) 6»(c - w 3 ) / du 3 f(u 3 , v 3 ) 

l-a-b 

n-l+a-b r b ~ 1± p- 
+ / du 3 6(u 3 + c)0(c- u 3 ) dv 3 f (u 3 , v 3 ) 

J-l-a+b J-b-l+pi 
r-l+a+b [b- 1 ^ 

+ / du 3 6(u 3 + c)e(c-uj 3 ) ' dv 3 f(cu 3 ,is 3 ). (A4) 

J-l+a-b J- a+ ±±^L 

Analogously for b > a the integral ( |Al| ) can be rewritten as 

-l+a-b ra+^p- 



I(a,b,c)= / du 3 9{u 3 + c) 9(c - u 3 ) / dv 3 f(u 3 ,v 3 

J-l-a-b J-b-l+p. 

+ / du 3 6(uj 3 + c)6(c-uj 3 ) dis 3 f(uj 3 ,v 3 ) 



-l-b+a 



-1+a+b rb-^pt 



+ / du 3 6(uj 3 + c)6(c-Lo 3 ) / du 3 f(u 3 ,u 3 ). (A5) 

-l+b-a J- a+ l+p. 



In order to be able to perform the integral ( |Al| ) we need to specify the intersection of 
three stripes ( |A3|) . The results ( |A4|) and ( |A5| ) are obtained for arbitrary positive a, b and 



c. Let us now take into the account that 

a = l-\Px\; 6 = 1-|^|; c=l-\(5 3 \ (A6) 

with \/3i\ < 1 and ft + ft + fa = 0. 

• Let us first consider the case when fts belong to the domain D\ ([0]). In this domain 
we have = \fa \ + \(3 3 \ and thus a = b + c — 1. So in the domain Di the following 
inequalities are respected: 

< a < 1 , < b < 1 , 0<a-6+l<l. (A7) 
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One may check that these inequalities result in 

c>-l + a + b; -c = -l + 6-a; and b>a. (A8) 
Thus employing (|A5|) we get 



-j — i, r b- 1 : p 



c 



J ( a Ae)| Dl = / rfw 3 / dv 3 f(u 3 , u 3 ). (A9) 



Analogously, in the domain D 2 we have \(3 2 \ = + \(3 3 \ and thus b = a + c — 1. 
The inequalities 

0<a<l, < 6 < 1 , 0<6-a + l<l (A10) 

result in 

c>-l + a + 6; -c = -l + a-6; and a>b. (All) 
Employing ( |A4| ) we get 

I(a,b,c)\ D = / du 3 / du 3 f(u 3 ,u 3 ). (A12) 



Finally, let us consider the case when belong to the domain D 3 . In this domain 
we have \(3 3 \ = + l/^l an d hence c = a + b — 1. Thus in the domain D 3 the 
following inequalities are respected: 

0<a<l, < 6 < 1 , 0<o + 6-l<l. (A13) 

One may check that in this domain 

-c> -l + b- a: -c>-l + a-b; c=-l + a + 6. (A14) 

Thus independently of a > b or a < b the integral over the intersection of three 
stripes ( |A4| ) or ( |A5| ) is again reduced to 

r-l+a+b fb-^ 

I(a,b,c)\ D = du 3 dis 3 f(uj 3j is 3 ). (A15) 

> I I 1+^3 
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B. CASE f < 



For completeness in this Appendix we present the result for ttN TDA H(w,v,^) in 
the ERBL-like and DGLAP-like type I and II domains for the case — 1 < £ < which is 



useful e.g. for NN — > 777* in the forward region |16| 



For w and v outside the domain — 1 < w < 1 and — 1 + |£ — £'| < t> < 1 — |£ — £'| 
the integral vanishes. 

For w G [-1; f] and u e [£'; 1 - f + f] (DGLAP-like type II domain): 

H{w, v, i ) = - / da / dp Fa, p, — A — £ . B16) 

For we[-l;(] and u G [— £'] (DGLAP-like type I domain): 

H(w, v,i) = - da dpF(a, p, — -— , — ^ . (B17) 

For we[-l; f] and v G [-1 + £' - & -£'] (DGLAP-like type II domain): 

H(w, v,0 = - / da / dpF(a, p, — — ^ . B18 

s J H ■'riq s s 

For tw G [£; -£] and v G [-£'; 1 - £' + £] (DGLAP-like type II domain) the result 
coincides with (|B1( 



w G [£; -£] and v G [£'; -£'] (ERBL-like domain): 

H(w, v,i) = - da / dp F a, p, ^— -, — ^ . (B19) 

io G [£; -f] and v G [-1 + £' - f; f] (DGLAP-like type II domain): the result 
coincides with (IB 181). 



w G [-£; 1] and v G [-£'; 1 - £ + £'] (DGLAP-like type I domain): 

/ i-{2 / 2^1-5 <?/; _ rr V — O 

H(w,v,0 = - 2 da dpF(a lP , ^— -, — I). (B20) 
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w e [-£ ; 1] and v e [£'; -f] (DGLAP-like type II domain): 

H(w,v,t) = -= da dpF{<T,p,^-^-,-A. (B21) 

w G [f ; 1] and v E [-1 + f - f; £'] (DGLAP-like type I domain): 

H(w, «, = H / <W , A ^> ^) • (B22) 
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